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Abstract. Elastic scattering of charm (c) and bottom (b) quarks via D- and B-meson resonance 
states in an expanding, strongly interacting quark-gluon plasma is investigated. Drag and diffusion 
coefficients are calculated from an effective model based on chiral symmetry and heavy-quark 
effective theory, and utilized in a relativistic Langevin simulation to obtain transverse-momentum 
spectra and elliptic flow (V2) of c- and Z?-quarks. The hadronization to D- and Z?-mesons is described 
by coalescence and fragmentation, and the resulting decay-electron spectra are compared to recent 
RHICdata. 

\o 

Introduction. Recent experimental results at the Relativistic Heavy-Ion Collider 
(RHIC) have given convincing evidence for the creation of dense partonic matter with 
large collectivity and opacity. A key challenge in the description of this strongly in- 
teracting quark-gluon plasma (sQGP) is the understanding of the microscopic reaction 
^ \ mechanisms, leading to its approximate behavior as a nearly perfect fluid. 

Heavy quarks (HQs) are valuable probes for the properties of the dense matter pro- 
duced in heavy-ion reactions, since they are expected to be created in the early stages 
of the collision. Recent measurements of the transverse-momentum (pj) spectra of non- 
photonic single electrons (e ± ) at RHIC, attributed to the decay of D- and 5-mesons, 
show a surprisingly small nuclear modification factor, [1, 2, 3], and large elliptic 



X 



flow, v e 2 [4, 5, 6]. To explain these findings, especially the large v 2 , quark-coalescence 
models [7, 8, 9] require that charm quarks are in approximate thermal equilibrium with 
light partons. A large degree of c-quark thermalization is, however, not supported by ap- 
proaches based on perturbative Quantum Chromodynamics (pQCD), e.g., using radia- 
tive energy-loss [10, 11]. While at lower pj elastic scattering processes parametrically 
dominate the energy loss (~ 1/ yfoQ) [12], a c-quark Raa compatible with the observed 
R% can only be obtained with unrealistically large values of the strong coupling con- 
stant [12]. Also the combined effects of elastic and radiative energy loss may not explain 
the experimental findings [13]. 

In this talk we introduce D- and 5-meson like resonance states in the sQGP [14] medi- 
ating elastic rescattering for heavy quarks. Employing pertinent drag and diffusion coef- 
ficients within a Fokker- Planck approach [17], we calculate HQ distributions in a flow- 
ing thermal QGP to simulate semi-central Au-Au collisions at RHIC [15]. Hadronization 
to D- and 5-mesons is described by a combined quark-coalescence and fragmentation 
model, and subsequent semileptonic decay electron spectra are compared to recent data. 
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FIGURE 1. Nuclear modification factor, Raa (left panel), and elliptic flow, V2 (right panel), of semilep- 
tonic D- and B-meson decay electrons in b=l fm, ^/snn = 200 GeV Au-Au collisions assuming different 
elastic HQ interactions in the QGP with subsequent coalescence, including the thermal weight factor 
described in the text, and fragmentation hadronization, compared to PHENIX and STAR data [1, 2, 5, 6]. 



Heavy-quark rescattering in the QGP. Lattice QCD (1QCD) computations of 
hadronic correlators and lQCD-based effective models suggest that mesonic reso- 
nance/bound states survive in the QGP up to temperatures of ~ 2T C in the light- and 
heavy-quark sector [18]. We here assume that the lowest pseudoscalar D- and 5-meson 
states persist above the heavy-light quark threshold [14]. Chiral and HQ symmetry 
imply the degeneracy with scalar, vector and axial- vector states. Pertinent resonant 
Q-q cross sections are supplemented with leading-order pQCD processes [16], using 
a s = g 1 1 (An) = 0.4. The evaluation of drag and diffusion coefficients within a Fokker- 
Planck model [17] results in HQ thermalization times which are lower by a factor ~ 3 
compared to pQCD scattering [15]. 

These coefficients are used in a relativistic Langevin simulation [12] for the rescat- 
tering of HQs in an isentropically expanding QGP fireball corresponding to b = 7 fm 
Au-Au collisions at RHIC. The expansion parameters are determined to resemble the 
time evolution of radial and elliptic flow in hydrodynamic models [19], with an ideal 
QGP equation of state with 2.5 flavors and a formation time of 1/3 fm/c (initial temper- 
ature Tq = 340 MeV). The proper thermal equilibrium limit in the Langevin process is 
implemented via the Hanggi-Klimontovich realization [20], with longitudinal diffusion 
coefficient B\ = TEA [12] (Einstein's dissipation-fluctuation relation). 

The initial HQ-p^-distributions and the relative magnitude of c- and £>-quark spectra 
are determined by fitting experimental D and D* spectra in d-Au collisions [21]. The 
corresponding e spectra saturate data from p-p and d-Au for p e T < 3.5 GeV [21, 22] 
with the missing yield at higher pj attributed to 5-meson decays, leading to a cross- 
section ratio of o b i/o C c = 4.9 ■ 10~ 3 and a crossing of D- and 5-decay electrons at 
Pt — 5 GeV. 

Hadronization and single-electron observables. The c- and ^-quark spectra from 
the Langevin simulation are used in the coalescence model of Ref. [7] with light-quark 
distributions from [23]. Here we take into account the thermal weight factor for the 
production of D* mesons relative to D mesons, (/no* 3 / 2 exp[— (mo* —mo)/T], 
as described in [7]. This leads to a reduced fraction of c-quarks which hadronize to 
D* mesons via coalescence, compared to our analysis in [15]. To conserve c- and b- 
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number unpaired HQs are hadronized via 5-function fragmentation. Finally, the single- 
e ± are obtained from D- and 5-meson three-body decays. Fig. 1 shows that resonance 
scattering leads to a substantial increase in v\ and decrease in R^, as compared to pQCD 
rescattering alone. Coalescence further amplifies v\ but also increases R^. The 5-meson 
contributions reduce the effects for pj > 3 GeV. 

Note that the nonperturbative resonance formation mechanism employed in this work 
importantly resides on a finite (equilibrium) abundance of (anti-) quarks, while perturba- 
tive calculations [11, 10] typically assume a maximum of color charges entirely residing 
in gluons. 

Conclusions. Assuming the survival of D- and 5-meson resonances in the sQGP, 
we have evaluated c- and Z?-quark spectra in an expanding fireball at RHIC within a 
relativistic Langevin simulation. The elastic resonance rescattering of c-quarks leads 
to an Raa down to 0.2 and V2 up to 10%, while ^-quarks are less affected. The HQs 
were hadronized in a combined quark-coalescence and fragmentation model followed 
by semileptonic D- and 5-meson decay. The resulting R^ and v\ are in reasonable 
agreement with recent RHIC data, suggesting that HQ-interactions via resonances may 
play an important role in the understanding of the microscopic properties of the sQGR 
especially the rapid thermalization of heavy quarks. 
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